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We describe general-relativistically a spherically symmetric stationary 

C>0 ' fluid accretion onto a black hole. Relativistic effects enhance mass accretion, 
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' in comparison to the Bondi model predictions, in the case when backreac- 
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tion is neglected. That enhancement depends on the adiabatic index and the 
asymptotic gas temperature and it can magnify accretion by one order in the 
ultrarelativistic regime. 

PACS: 04.40.-b, 98.62.Mw, 04.20.-q 

I. INTRODUCTION. 

In this paper we re-examine the spherical gas accretion onto a black hole, paralleling 
previous studies of fluid accretion of Michel and Shapiro- Teukolsky ( 0). It is shown 
that relativistic effects can lead to a bigger mass accretion than that predicted by the 
corresponding Bondi model 

The order of this paper is as follows. Next section presents spherically symmetric Einstein 
equations expressed in the language of extrinsic curvatures. A suitable choice of a gauge 
condition leads to a "comoving coordinates" formulation that is particularly suitable 
for the description of selfgravitating fluids. In Section III we show that the original set of 
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integro-differential equations can be reduced to an integro-algebraic problem, whose solution 
would constitute a new stationary, general-relativistic solution of selfgravitating polytropic 
fluids. That model is complete - it includes the back effect exerted by matter onto a metric 
- therefore it is capable to describe a stationary phase of the interaction of (even) heavy 
clouds of gas with a relatively hght center. Section IV discusses a case when backreaction 
can be neglected. Under some circumstances, an accretion is described by a set of purely 
algebraic equations. Section V proves several quantitative and qualitative properties of 
accreting solutions. It is shown that the Bondi model relation between the asymptotic 
and sonic speeds of sound appears as a limiting case of relativistic formulae. Section VI 
compares predictions of the Bondi model and of the relativitic solution without backreaction. 
Relativistic magnification of the mass accretion becomes noticeable in the case of infall of a 
hot gas, when the correction factor can be bigger than 2.4(1 + ^), where T is the polytropic 
index. 



II. EQUATIONS. 

We will use a spherically symmetric line element 

ds^ = -N^dt^ + adr^ + R^dO^ + sin^ (2.1) 

where N, a and R depend on t (asymptotic time variable) and a coordinate radius r. we 
will work in extrinsic curvature variables. Thus we need the mean curvature of centered 
two-spheres in a Cauchy slice, 

and the extrinsic curvature 

dtiV^R^) 



trK 
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N^R 

Kt-K^e-^-litrK-K). (2.3) 
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Let be the energy- momentum tensor of matter fields, p = — Tq and > = NT^, R^y 
be the Ricci tensor and R the Ricci scalar. 

The Einstein constraint equations Rq^ — gQ^R/2 = SttTq^ can be integrated to yield, 
assuming asymptotic flatness. 



Rp = 2Jl- — + —j^ R^pdR + T, (2.4) 



Ci-^^k^'-^dR XtrKiedR 



RKl - RtrK = ^^^2 ' (2-5) 

where m is the asymptotic mass and 

Impose the integral gauge condition 

, ^ i^mtEzJSlf (2.7) 

where r is given in ( p.6|) . One can show that in vacuum ( p.7|) is satisfied identically. 
Differentiation of (12^71 ) with respect r yields, after some algebra, 

R{tTK-Kl)^^^^ = Q (2.8) 

which implies 

JV = = (2.9) 

in geometries without minimal surfaces and with trK ^ K^. Thus in this gauge coordinates 
are " comoving" - each particle of matter carries a fixed value of a radial coordinate " r" . 
The energy-momentum tensor of a self gravitating fluid reads, in comoving coordinates, 

T^u = {p + P)U^U,+Pg^,. (2.10) 

Here U^JJ^ = —1. Notice that the presure is p = T!^ = Tg. This space-time foliation is 
regular even at the vicinity of the boundary of a black hole, in contrast with other approaches 
im, P in which the Schwarzschild geometry is foliated by polar gauge slices. 



Define a mass function 

miRir)) = 27r f R^ppdx (2.11) 

^0 



'0 

wliere R is an areal radius. Tlie mass evolves as follows 



dom{R{r)) = -2Ti[NR^trK - K;)p]{r) (2.12) 



Direct differentiation of m{R) gives 

drm{R) 



= 2rfpp. (2.13) 



The remaining relevant equations are the two continuity equations 



Ndrp + drN{p + p)^0, (2.14) 



dop^ -NtrK{p + p). (2.15) 



and the Einstein evolution equation 
dtiK - trK) = ^{Kf - ^ - -^drN + — + SttNT; + ^N{trKf - ^trKK; 

(2.16) 

The rate of accretion m of mass along orbits of a constant areal radius R is equal to 

m{R) = {do - RdR)m{R(r)) = -AnNR^Uip + p). (2.17) 

where 

U = doR/N^^{trK-K;). (2.18) 

III. STATIONARY DESCRIPTION OF A SELFGRAVITATING FLUID. 

All results of this section hold true for systems with collapsing or exploding matter. 
Assume a compact cloud of a fluid. We will say that an accretion (explosion) is stationary 
if 



i) the mass accretion 

m = {do- RdR)m{R{r))\R=const 

on a central body is constant in time; 

ii) the radial fluid velocity U = is constant at a fixed value of the areal radius i?, 
(9o - Rdn)U = 0; 

iii) the energy density at a fixed areal radius does not change in time; 

iv) asymptotically, i. e., close to the outer boundary of the collapsing (exploding) gas, 
its speed U is much smaller than the speed of sound = dpp, Uoo ~ 0. (Expanding fiuid, 
in turn, would be subject to a condition U ^ at the inner boundary 

At first, we shall prove the following fact 

Theorem 1. Under conditions i-iii), rh does not depend on R within the fiuid filled 
zone, dRrh = 0. 
Proof. 

Equations (g^) and (^181) yield 

dniUR^) = RHrK. (3.1) 
/^From rri{R) = —AnNUR'^lp + p) we obtain durn = I -\- II + III, where 

I = N{p + p)dR{UR^) 

II = NUR^Orp + UR^p + p)dRN 

III = NUR^dnp. (3.2) 

Using ( |3.1| ) one writes / = N R^trK{p + p), while the stationarity condition ii) allows 
one to write /// = R^d^p = —N R^trK{p + p) (the second equality follows from (|2.15|) ). 
Thus / + /// = 0; since 11 = (due to the momentum conservation ( p.l4|) ), we arrive at 
dRm{R) = 0. 

Assume the equation of state 

p = Kp'^, (3.3) 
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r being a constant and define the speed of sound as = dpp. We assume tliat 1 < T < 5/3, 
since we are primarily interested in comparing predictions with the Bondi model, but it is 
quite likely that much of the forthcoming analysis applies to adiabatic indices in the standard 
in astrophysics range (1, 2). 

Let us point out that astrophysicists use a different equation of state, p = CvF (where 
n is the baryon number density); that reads in our notation 

p = C X exp(r j dpj^j^). (3.4) 

Both approaches agree for F = 1 or in the newtonian limit when F ^ 5/3, but they disagree 



for F = 5/3. The momentum conservation equation (|2.14|) can be integrated, 



= + (3-5) 

where k = and the integration constant is equal to the asymptotic speed of sound 
of a fluid. 

( |?75| ) asymptotically {m/R « 1) yields the Bernoulli equation, hence it can be regarded 
as the general-relativistic version of the latter. 

/^From the relation between pressure and energy density, one obtains, using equation 

m 

P = Poo[alaoo) =Poo[—^+ V-', (3.6) 

where the constant poo is equal to the asymptotic mass density of a fluid. /^From the evolution 
equation one gets, using the stationarity assumption, 

RdnU = h^pRfdnN - m{R)^ - AnRNp. (3.7) 



Equations ( |2.4| ) and ( |2.5| ) give 



UdjiU = ^dnipR) - ^ + ArcpR. (3.8) 



Comparison of (|3.8|) with ( p.7|) yields an ordinary differential equation 



(p+py, 



integration of that, with the asymptotic condition at spatial infinity 
the following relation between the lapse N and the mean curvature pR 

pR 



pR 
2 



(3.9) 

1 leads to 

(3.10) 



where 



(3{r) = e 



lQ7T(—p—p)—h-ds 

J r s 



(3.11) 



The substitution of trK (as calculated from the continuity equation ( p.lSj) ) into (|3.1|) gives, 
employing the stationarity condition, 

dRP 



OrHIUIR' 



p + p 



(3.12) 



dnip+p) _|_ ag. 



The last term can be presented in another form (due 



Notice that -- ,^ , _^ . 

p+p p+p p+p 

to relations between a^,p and p), |^ = ^^9/? ln(a^/r + 1). That leads to the following 
solution of ( |3.12| ) 

(aVr + 1)1/(^-1) 



u = c- 



(3.13) 



i?2a2/(r-i) • 

The whole set of equations describing the collapsing stationary fiuid is given by ( |3.13|) and 
the previously written equations (|3.5| , p.lO|) and (|3.11| ). Calculation of dR\n.{o? + F), with 
given by (|3.5|) and being specified above, yields 



(3.14) 



One easily obtains from (|3.13|) that 

1 



dniU'R'' 



(3.15) 



The insertion of (|3.15|) into ( p.l4| ) gives 



dR{U'R%l 



_ 16R m{R) 
a2p2ji2j - 2R 
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+ 2-nR^p - U 



(3.16) 



We define sonic points as such where the the equahty U = holds true. Let be a 
radius of a sonic point. Equation ( p.l6|) yields the relation 

al{l-^^ + c.) = U! = ^ + c., (3.17) 

where = a1 = a^i^R^), m^, = m{R^) and = U'^{R^). 

The constant C in formula ( |3.13[ ) can be expressed in terms of a*, [/* and i?*, that is as 
a function of c*, m(_R*) and R^. The infall velocity U reads 



t^ = f^*7i(f^ (3-18) 



Above t/* means a negative square root in the case of falloff towards a gravity centre and a 
positive square root in the case of exploding gas. 

The rate of accretion of mass ( p.l7| ) can be conveniently expressed by characteristics of 
the sonic point i?*, 

9 / 3m^, \i/2 /a(i?*) \ 2/(r-i) ^ a^^^, ^ 
m = -47ri?2p^(l-— ^ + c,) f/,(^^) + (3.19) 

For the sake of completeness we write down the space-time line element with the areal radius 
chosen as the radial coordinate, 

ds^ = dt\-N^ + —4^) - Af3—dtdR + —\^dR^ + R^dn\ (3.20) 

[pRy pR [pRy 



IV. RELATIVISTIC ACCRETION: NEGLECTING BACKREACTION. 

The quasi-stationary accretion shall apply to the description of black holes interacting 
with a fiuid. The description of the accretion onto other compact bodies (say, neutron stars) 
is more complex, since there can appear shocks, that are excluded in our picture. The above 
model can be valid only if shocks are absent, for instance when the inner boundary of a 
collapsing shell of gas is disconnected from the surface of a compact body. 

All hitherto proven results are exact and - under the preceding reservation - they refer to 
a fully nonlinear stationary system consisting of a central mass and a cloud of gas that would 
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dynamically influence a geometry through a backreaction. If the gas is heavy, comparing 
with the central mass, then l3{R) is nonconstant; metric functions do depend on the infalling 
matter. That means that backreaction should be taken into account in description of such 
a system. 

If the contribution of a fluid to the total asymptotic mass of a system is negligible, i. e.. 



TTlf 

and p < p, then P ^ 1 and 



/ drr'^p « m (4-1) 

jR>2m 



iV«f =«f -l^ + C/'^. (4.2) 

That would suggest that in this case the standard schwarzschildean metric constitutes a 
valid approximation. There is, however, one subtle point. The reasoning of the former 
section shows that in order to neglect the effect of backreaction the following condition 

c, =2rfp<<^. (4.3) 

must hold at a sonic point. That can be interpreted as the demand that not only is close 
to pR/2 but also OrN shall be approximated by dR{pR)/2. 

We will say that backreaction is negligible if both conditions (|4.1| ) and ([4.3|) hold true. 
In such a case m ~ m^, and one obtains at the sonic point 

alii - 1^) = = (4.4) 

The remaining two equations describing accretion are 

Rim/ 1 N2/(r-i rN2/(r-i) 



I 1 N2/(r-i / r N2/(r-i) 
and 

r + 

= -r + ("-6) 

This is a purely algebraic system of equations, describing the fluid accretion in a flxed 
space-time (Schwarzschild) geometry. 
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V. RELATIVISTIC ACCRETION WITHOUT BACKREACTION. 



Numerical analysis demonstrates - as pointed first by Michel [|l| - the existence of two 
branches of solutions of the relativistic fluid equations. An analytic proof is given below. 

In the first part we prove the existence of a sonic point in a black hole spacetime endowed 
with a Schwarzschild metric. That black hole - fluid system is shown to possess a sonic 
point; that leads, through a construction outlayed in the second step, to the existence of 
two accreting solutions. 



A. Sonic points 



Define 



L = a' + r 

+r 

^ r . n(r-i)/{2r) y'^--^) 



where is given by ( [4.5D with parameters and specified by 

The equation L{R^) = P{R^:) for a sonic point can be written as 1 + y{3T — 1) = 
'^i'^oo + r)?/('"+^)/(^'"), where y = 1 — 3m/ (2i?*). One has to demand that y > (i. e., 
i?* > 3m/2), since at y = (or R^, = 3m/2) the coordinate system breaks down. Notice a 
numerical mistake in [|l| which led Michel to the wrong claim that sonic points must exist 
outside a sphere of a radius 6m. In fact they may exist even inside a black hole, although - 
as we point out below - that would contradict established views on properties of matter. 

The left hand of the equation in question is bigger than its right hand side at y = while 
at y = 1 the opposite holds true. Since both sides are continuous in x, their graphs must 
intersect somewhere. Since 1 + ?/(3r — 1) increases at a lower rate than 3(a^ + T)y^^^^^^^'^^'' 
for r < 5/3, there exists a unique sonic point characterised by y^:. The case with y^ < 1/2 
(i. e., when R^: < 2m) is physically noninteresting. In that case the speed of sound would be 
bigger than the velocity of light and the dominant energy condition [^] would be broken, even 
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outside of a black hole. One easily infers that y* is a monotonously decreasing function of the 
asymptotic sound density a^. Therefore there exists a critical value of which separates 
solutions that are subluminal from unphysical solutions that become superluminal. 

An interesting feature of the Bondi model is the simple relation a^/a^ = 2/(5 — SF) 
for r < 5/3. Below we will show that this relation appears in the nonrelativistic limit of a 
relativistic formula. 

Theorem 2. Let and al be the asymptotic and sonic speeds of sound, respectively. 
Define q;= (r-l)/(2r) 

^ = a((9r-5)ln4-1.5) 

and 

B ^^a'^T{9T -7). 

If sonic points are exterior to a black hole then 

1 a? 1 , , 

5-3r I p aj — fj2 — 5-3r , a aj ^ ' ' 

Proof. 

Define x = ^ (or alternatively, x = a^/(l + 3a^)), and 

5 — sr 

F = x{l- 2>xy + r(l - 3a;)i+" - r(l - 3x) ^x - Bx^ 

5 _ 3P 

* = x(l - ?,xY + r(l - 3x)^+°^ - r(l - 3x) — X - Ax'^. (5.3) 



The sonic point equation can be written as 

al 2 X 

Ax+-. (5.4) 

X 

It suffices to show that F > and ^ < 0. We shall deal with the first inequality. The 
second derivative of F with respect x reads 

F"^Z^^{l-?,xr-'G{x), (5.5) 
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where 

QF — 7 V -\- 1 1 

G{x) = '-^{1 - 3x) - 3x^ - -(r - l)(9r - 7)(1 - 3x)^-". (5.6) 

One shows that G' < [ST — 5)/r < 0; thus G{x) is decreasing for < x < 1/4 and 
1 < r < 5/3. Therefore if F"{xq) = then F"{x) < for any x > Xq. That means, taking 
into account that F"{0) > and F'{0) = 0, that if F' vanishes at a point xi, then it must be 
negative in the interval (xi, 1/4). In conclusion, either F is increasing (and then it achieves 
its minimum at x = 0) or it has a single extremum (a maximum) in (0, 1/4). Notice now 
that F{0) = 0. Thence in order to show that F{x) > it is enough to show that F{x) is 
nonnegative at x = 1/4, when 

nm = (r + i)^ + ^- ^- (r - i)^(9r - 7). (5.7) 

A numerical calculation shows that -F(l/4) > and the equality is achieved only at F = 1. 

In a similar vein one proves the other inequality \E' < 0. At x = one has \E' = 0. On 
the other hand, 

QF — 5 / 3n 

^' = (1 - (1 - 3x)")-3xa(l - 3x)"-^ - 2ax((9F - 5) ln4 - ^). (5.8) 

Employing the estimate 

1 - (1 - 3x)" < 4axln4 (5.9) 

which is valid for < x < 1/4 and 0.2 > a >0, one arrives at < 0. Thus the function 
is nonnegative, as desired. That ends the proof. 

Let us remark that (|5.2|) can be written as, resolving the biquadratic inequalities. 



60 - J6i + 4a^52 



26o 



>al> 



60 - J 51 + Aa'LSi 



(5.10) 



where 
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(5.11) 



Asymptotically, i. e., for x one obtains the Bondi equality al/a^ = 2/(5 — SF) for 
r < 5/3. If r = 5/3 then the above gives asymptotically 1.12 x Ooo > > O.Saoo, in a good 



agreement with the exact formula al = y 5/6aoo- If a sonic point is located at a horizon of a 

black hole (that is, a* = 1) then ( p.2|) (or the above inequalities) yields 0.79 > aoo > 0.5 for 
r = 5/3. Notice also a rough bound al > 1.6 x which is valid for any T and a^; outside 
of a black hole al < 1, therefore one infers that the asymptotic speed of sound is less than 
1. 

Let us point also that equation ( ^.4|) implies that asymptotic sonic points can exist only 
for models with adiabatic indices T < 5/3. 



We show that at least two solutions {a{R),U{R)) bifurcate from i?*. Define (2q- as a 
solution of the equation 




B. Existence proof 





^From that and (U) follows that f/^ = U^^JrJr and 

2 T{R,/R)f^ 



(5.13) 



5-iRjR)f'' 



where /? = |(r - 1) and 5 = 1 + T/al 



A straightforward calculation gives 



_d_ 

dR 



In L{aa) 



P{r./rY 



(5.14) 



R{5 - {R,/RY) 



and 
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d , , 2(r- l)i?,f/2(l 



2n V R/R* 



— In P(a,) = ^ « ' . (5.15) 

dR ri?2(l + t/2(3 - 4f + )) 



L and P are equal at i? = i?* and they are decreasing in the vicinity of R = R^,. Morever, 
djiL = OrP at the sonic point i?*. A careful investigation shows, however, that second 
derivatives are both locally positive and 

rf2 £ 29/14 + 7a2/2 

\nPM\R=R, = -—\nLM\R=R^ X ^' ^ ^ * ■ (5.16) 



One observes that P{aa)\R=R, > -j^ L{aa)\R=R, if F < 79/49. This reasoning can 

be valid for adiabatic indices F < 5/3 assuming that the exponent 7/2 in ( |5.12| ) is replaced 
by a;e(-oo, 4.5 - v^) 0- Therefore OrL < OrP, for R > R^ and OrL > OrP for R < R^. 
Thus locally P > L. 

On the other hand, notice that L(a^ = 0) > P{a^ = 0) and L(a^ = oo) > P(a^ = oo), 
for all values of R. 

L and P are differentiable functions of their arguments. Combining the above facts one 
infers that, due to the continuity of L and P, there must exist at least two solutions in a neigh- 
bourhood of -R*. Those solutions coincide aX R = R^, due to the above construction. The set 
of those points constitutes at least two branches. Since da?{L — P) = 1 — AU"^ / (jP'R^a?) ^ 
at any point of a solution branch with R R^,, the implicit function argument would be 
used to extend the interval of the existence onto a whole bounded domain. Those solutions 
are differentiable for R ^ R^... 

One of the solutions is supersonic below -R* and subsonic above i?* and it can be inter- 
preted as describing collapse of matter onto a black hole. The other solution is subsonic for 
R < R^: and supersonic above; it can correspond to an exploding gas. 



C. Qualitative results 

In what follows we shall deal with a solution that is subsonic asymptotically, i. e., 
describes accretion of a fluid. 
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Theorem 3. An asymptotically subsonic solution of the system ( |4.2| - [4.6|) satisfies 
following conditions: 

i) If i? 7^ R^: then Or^U'^R^) > and the speed of sound decreases, dno^ < 0, with the 
equality only at spatial infinity; 

ii) > ^ for < R, and < ^ hi R > R,. 

iii) Inside the supersonic region a?{pRY /A < < 2m/R. 

iv) Mass density p monotonously decreases and p is bounded in the supersonic region, 
R < -R*, 

4:771 1 

P<Poo[i + (r-i)-^]— . (5.17) 

The proof is postponed to the Appendix. 

The estimates of ii) and iii) in Theorem 3 require an explanation. It proved to be 
convenient to define a sonic point by requiring that a^(pi?)^/4 = f/^ instead of the condition 
(used in the Bondi model) o? = U'^. Therefore the speed of sound can be bigger than infall 
velocity in regions close to horizons if the factor pR/2 is significantly smaller than 1. In 
the traditional terminology such a solution would be called subsonic. Numerical data of the 
next Section show that the value \U\/a at a horizon depends strongly on the location of a 
sonic point - on the ratio R^/{2m), which in turn depends on the asymptotic speed of sound 
Ooo- \U\/a decreases with the increase of the asymptotic speed of sound. 

VI. BONDI MODEL AND THE RELATIVISTIC SOLUTION 



The insertion of ( [4.4|) and (|4.6| ) into (|3.19|) (with = 0) yields the mass accretion rate 



777 = irm 



oPoo / a? \ (5-3r)/2{r-i) ^ a^x. 9x , 

'{l + f){l + 3al). (6.1) 

One can write that as 

777 = QmB, (6-2) 

where 
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2 Poo/ 2 N(5-3r)/2(r-l) 

^B = ^m-[^-^) (6.3) 

is the mass accretion rate predicted by the Bondi model and 
^ /(5 - 3r)a2x (5-3r)/2(r-i) 

^ = ( .2 ) (l + 3a,^)(l + -^). (6.4) 



VL can be interpreted as the relativistic correction factor. 
Apphcation of Theorem 2 leads to useful estimates for Vt. 
Theorem 4. Assume 1 < F < 5/3. The relativistic correction factor satisfies 

(l + 3a^)(l + ^) 
>VL> 

(l + 3a^)(l + ^)e-^, (6.5) 



where 



Proof. The definition of Vl and Theorem 2 yield immediately the forthcoming inequali- 
ties, 

(l + 3a,2)(l + 



I) 

> 



2B 3ai \(5-3r)/2(r-l) 



^ + (5-3r) 

n > 

(1 + ?{£zi)c)(5-3r)/2(r-i)- y'^-') 

2 

Bounding from above the left hand side of ( |6l7|) by (1 + 3a^)(l + ^) yields the first bound 
of Theorem 4. The proof of the other inequality bases on the obvious estimate 

(1 + 2(r-l) ^^(5-3r)/2(r-i) ^ (g_8) 



Taking into account Theorem 2 and its implications stated in ( [5. 101 ) and below formula 
( ^.11| ), one can write the above in terms of asymptotic data, 
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1 + 3- 



6n- 



52 + 4a|^^^^_^5o_ 



26o 



(6.9) 



The relativistic correction factor Q is close to 1 when << 1, i. e., when the asymptotic 
gas temperature is low. Q is bounded from below by 0.99. (|6.5|) yields, in the ultrarelativistic 



regime 



4(1 + ^)>^^> 1.6(1 + 1; 



(6.10) 



Ultrarelativistic effects enhance accretion, with the strongest effect for the isothermal gas 
with r = 1. The enhancement is smaller for F ^ 5/3, as seen from the preceding estimate. 

The Bondi model fails only in describing the hot gas mode. The correction factor Q tends 
quickly to 1 when sonic points are far away from the Schwarzschild sphere. For instance, if 
F = 4/3 then < 7 at Rjm = 2 but f2 < 1.1 at Rjm = 25. 

We analyse numerically a relativistic gas, with the adiabatic index F = 4/3, falling onto 
a black hole. Results complement analytic estimates and they are comprised in the table. 



R,/i2m) 


a'iR.)/al 


* aiRn) 


* U{Rn)) 


*p{Rn))/poo 


P{R*)/Poo 


p{2 X R^)/poo 


500000 


2 


0.168 


0.999 


1.41x10^ 


8 


3.14 


50 


1.99 


0.173 


0.92 


1595 


7.89 


3.9 


5 


1.91 


0.34 


0.78 


56.4 


6.95 


3.46 


1.1 


1.62 


0.86 


0.53 


4.82 


4.24 


2.1 



* Rn = 1.001 X 2m 



Some features of accreting solutions depend in a crucial way on the location of sonic 
points. When sonic points are close to a horizon, the speed of sound is close to one while 
the infall velocity at a horizon is smaller and it barely exceeds 1/2. When sonic points are 
far away from a horizon, R^, » 2m, the infall velocity nears to the speed of free fall (f/ ~ 1 
close to a horizon) while the speed of sound is then much smaller than U. An interesting 
fact is that the energy density changes quite moderately - by a factor of the order of unity 
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- if sonic points are close to the Schwarzschild sphere. In contrast with that, solutions with 
» 2m are characterized by a rapid growth - up to ten orders - of the energy density 
near the horizon. The energy density changes by a factor not greater than 8 in the region 
exterior to a sonic point with i?e(it!*,oo); that type of moderate decay is common for all 
solutions, irrespective of the value of 2m/ R^. That actually follows from Theorem 2, which 
bounds al/al^ and - consequently - also p{R*)/poo- Solutions with sonic points close to a 
horizon have a approaching 1 and they describe a high temperature (circa 10^° K) gas, with 
^oc ^ 0.5. 
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Appendix 



Proof of Theorem 3. 



1) Equation ( p.l6| ) yields, ignoring the backreaction term, the crucial relation 



dRiU'R'){l 



) 



IQR , m 




(6.11) 



p2 ^2R 



2) The first observation, that signs of d^a? and of dji{U^R) are opposite and that they 
vanish simultaneously at finite values of R, can be drawn from (|3.15| ). 

3) Let R^: be a position of the sonic point; thus f/^(-R*) = m/2R^:. Assume that in the 
vicinity of i?* the expression dR{U^R^) is strictly negative. Then ( |6.11| ) yields U^R < 171/2 
for R < R^: and U^R > m/2 for R > R^. Therefore U^R is increasing in the region of 
interest and that is incompatible with the assumption that dfi{U'^R'^) is strictly negative. 
Thence it must be weakly positive at least around the outermost sonic point. That in turn 
implies that in a neighbourhood of a sonic point 2U^R < m for R> R^ and 2U^R > m for 
R < R^, 

4) The expression d^lU'^R'^) cannot have zeroes. For, let it vanish at some Ri > R^,; 
( |6.11| ) gives U'^{Ri) = and we would have dfi{U'^R) > at Ri. But that is incompatible 
with the assumption that dniU'^R'^) = at -Ri. 

Let us now consider a region R < R^.. If _Ri is a zero point of dRiJJ'^R^) but the latter 
does not change sign at then 2U^R decreases for R < Ri towards the value m and 
increases for R > Ri (due to estimates proven in the final part of 3)), Hence dptiU^R) = 
at -Ri. But that contradicts OrIU^R'^) = at -Ri. Similarly, if dRllf^R^) changes sign in the 
vicinity of then we are led to contradiction. 

Thus dii{U'^R'^) > in the domain of existence of the solution. That implies, with 
conjunction with ( |6.11| ), that in the supersonic zone f/^ > m/{2R) and that f/^ < m/{2R) 
in the subsonic zone {R > R^:). That accomplishes the proof of ii) 

5) Rewrite Eq. ( |3.14D , with backreaction terms being dropped out. 




(6.12) 
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Let i? be a largest point R < such that U'^R = 2m; then dR{U'^R)\R < and from (|6.12|) 
follows 9i^a^|* > 0, in contradiction with hitherto proven monotonic falloff of a^. That 
shows the bounds of iii) - that U^R < 2m. ( |6.11[ ) and ii) imply, in the supersonic region, 
a^p^R^ /4: < 1. Extremal values of the speed of sound (achieved at a horizon of a black hole) 
cannot exceed A/{pRY while the speed of infalling particles does not exceed 1, the speed of 
light. 

6) The decrease of the speed of sound together with (|3.6|) lead to the conclusion that 
the mass density also decreases, Oro? < and Orp < 0. The numerical estimates of iv) are 
obtained from inserting inequalities proven in 5) to the expression ( |3.6|) . 

That ends the proof of Theorem 3. 
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